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Estimation of Faults in DC Electrical Power System 


Dimitry Gorinevsky*^, Stephen Boych, and Scott Polk 


Abstract — This paper demonstrates a novel optimization- 
based approach to estimating fault states in a DC power 
system. The model includes faults changing the circuit topology 
along with sensor faults. Our approach can be considered as 
a relaxation of the mixed estimation problem. We develop 
a linear model of the circuit and pose a convex problem 
for estimating the faults and other hidden states. A sparse 
fault vector solution is computed by using l\ regularization. 
The solution is computed reliably and efficiently, and gives 
accurate diagnostics on the faults. We demonstrate a real-time 
implementation of the approach for an instrumented electrical 
power system testbed at NASA. Accurate estimates of multiple 
faults are computed in milliseconds on a PC. The approach 
performs well despite unmodeled transients and other modeling 
uncertainties present in the system. 

I. Problem 

We consider a DC (direct current) electric circuit with 
sources, loads, and switching elements. Voltage and current 
measurements are available at certain circuit locations. The 
problem is to estimate, from the measurements, the fault 
states of the circuit. The faults are defined as the deviations 
from the nominal state. The source voltages or loads can 
differ from their nominal values; voltage and current sensors 
can be faulty; assumed open/closed states of relays and 
breakers might differ from their actual states; short-circuit 
and open-circuit faults are possible, including shorts to the 
ground. We assume that several faults might be present si- 
multaneously, however, it is known that most of the potential 
faults are not present. 

Let y £ be a vector of observations (measurements). 
The fault states are described by a vector f £ - : R Nf . The 
problem is: given the observations y, estimate faults /. The 
states of the circuit (e.g., voltages and currents) are described 
by a vector x £ $t Ns and have to be estimated along with /. 
Section III below describes a circuit model of the form 

0 = Ax + Bf + Z, (1) 

y = Cx + Df + rj, (2) 

where matrices A, B, C, D have appropriate dimensions. 
The process noise vector £ £ $t Na describes the modeling 
error; the measurement noise rj £ $R Nm describes the obser- 
vation error. State equation (1) and observation equation (2) 
are similar to linear state space model commonly used in 
estimation, except the model (1), (2) is static. 
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We estimate the faults / and the states x by solving a 
quadratic programming (QP) problem of the form 

minimize * || Ax + Bf\\ 2 Q + * \\Cx + Df - y\\ 2 R + X T \f\, (3) 

where ||z||q = z T Qz; Q, R are positive definite matrices; 

| /| is the component- wise absolute value; A is a vector with 
positive components. The last term provides a weighted l\ 
penalty for the unknown components of vector /. 

Optimization-based estimation of electrical power systems 
is well established area, e.g., see [1], [19] for an overview, 
more discussion is in Section II. The novelty of the formula- 
tion (3) is that the faults / are introduced in a structured way 
separately from the states x. Quadratic penalty is used for 
mismatches £, tj in model (1), (2), which are usually small; an 
1 1 penalty is used for the fault vector /, which is usually large 
and sparse. We demonstrate that such formulation allows 
for efficient and accurate estimation of the faults. Work on 
relaxed mixed estimation related to the proposed approach 
is discussed in [29], [30]. 

Problem (3) can be interpreted as Bayesian Maximum 
A posteriori Probability (MAP) estimation of x and / 
given y. The MAP formulation assumes Gaussian noises 
£ ~ AT(0, Q” 1 ) and y ~ N( 0, f? -1 ) in (1), (2) yileding 
the quadratic terms in (3). A Laplacian prior distribution 
assumed for the fault vector / yields the l\ term. No prior 
is assumed for the state x. 

The main benefit of using l-\ regularization in (3) is that 
for properly selected weights A the solution vector / is 
sparse. See [6], [9], [27], [28] for the recent work and 
further references on sparse solutions using l\ regularization. 
Another advantage of QP formulation (3) is the ease of 
including additional linear constraints. Section V includes 
positivity constraints for some voltages and currents. QP 
problem (3) can be solved in milliseconds for hundreds of 
states, using modern interior-point methods, see [5], 

Today electric utilities use EMS (Energy Management 
Systems) to monitor, control, and optimize the transmission 
and generation facilities. Optimization-based estimation is 
presently used for power systems monitoring, see [19], [1]. 
On-line QP optimization is also used in Model Predictive 
Control systems broadly employed in industry, see [24]. 

II. Motivation 

There is a large body of prior work in diagnostics of 
electric power systems. Much work is focused on specific 
electric power system units and elements, such as electric 
machines, motors, generators, inverters, batteries, solar cells, 
relays, and other. Work on diagnostics of power distribution 



systems is primarily focused on AC (alternating current) 
systems, e.g., [8], 

A few papers consider large vehicle electric systems: inte- 
grated diagnostics of international space station is discussed 
in [11]; diagnostics of marine vehicle power system in [16]; 
integrated diagnostics and prognostics of aircraft electric 
system in [12]. In most of this work, Al-type computational 
reasoning is used. 

Our formulation is related to the approaches in 
optimization-based estimation of power system state that 
were first introduced in 1969, see [4], [2] and the books 

[1] , [19] for an overview. The state vector is determined 
by least squares fit of the nonlinear equations relating the 
state variables (bus voltage magnitudes and phase angles) 
and the measurements (voltage and current magnitudes). The 
problems are generally nonconvex and an efficient solution 
with global convergence is difficult to achieve. 

Some of the power system estimation work explicitly 
considers outliers in the measurements data (an equivalent of 
the sensor fault) and unknown states of the circuit breakers 
(might be caused by a fault). Using l\ model fit error 
instead of the quadratic error for countering and detecting 
bad measurement data is discussed in [10], [14], also see 

[2] , [4]. The papers [25], [13], [7] look into the problem of 
determining network topology errors. In [25], [13], mixed 
integer problems for determining topology are formulated. 
The cited papers pursue nonconvex optimization problems. 

Using GPS for time synchronization made possible accu- 
rate phase measurements and harmonic state estimation of 
power systems, e.g., see [17]. Harmonic models are linear, 
which makes them close to the models in this paper. An 
approach using l \ penalty in harmonic source estimation is 
discussed in [15], 

This paper considers a DC power system with linear 
constitutive equations (1), (2). We show how the fault pa- 
rameters can be chosen to preserve the linearity, including 
the topology changing faults. Convex optimization problem 

[3] can be solved using standard QP software. Our approach 
can be extended to fault estimation in harmonic power system 
analysis, where the equations are linear. 

We implement the approach for the ADAPT EPS testbed. 
This testbed was developed and is maintained at NASA 
ARC as an experimental platform for research in integrated 
systems health management, diagnostics, and prognostics. 
More detail on ADAPT can be found in [23]. Earlier work 
in ADAPT diagnostics was based on Al-type reasoning and 
other heuristics, see [3], [18], The proposed optimization- 
based estimation approach works very well for ADAPT. 
Multiple faults are diagosed correctly in milliseconds. 

III. Model 

This section formulates a linear model of the form (1), (2) 
for a DC circuit. 

A. State equations 

We formulate the Sparse Tableau Analysis (STA) equa- 
tions of the nodal DC circuit in a standard way, e.g., see 


[22], STA models are the core of most standard small-signal 
linear analysis of electric circuits. This section extends the 
STA to include a linear model of the faults. We consider 
several types of faults. Some of them or all of them might 
be present and need to be considered. 

The circuit has N]\ r nodes; each node has a voltage e/, : with 
respect to the ground. The circuit contains Nb branches, 
each has current j) and voltage drop v/ . The signs of the 
currents and the voltage drops are relative to the directions 
of respective branches (graph edges). The incidence matrix 
G £ s ft N N,N B h as en t r j es g kl — i if thg branch l leaves 
the node k, g ki = —1 if it enters node k, and g^i = 0 
otherwise. The STA equations comprise Kirchhoff’s Current 
Law (KCL), Kirchhoff’s Voltage Law (KVL), and Branch 
Constitutive Equations (BCE), which are discussed below. 

The KCL equations for the currents j £ $R Nb are 

S B Gj = 0, (4) 

where Sb £ $R Nl ' Njv is a selection matrix that un-selects the 
boundary nodes of the circuit where KCL does not hold. 

The KVL relates voltage drops v £ ‘ft Y,i to the node 
voltages e £ 5ft Yjv 

G t e = v (5) 

Finally, the BCE relate the voltage drop and the current. 

-K I j + K v v = f G , (6) 

where Ky and Kj are diagonal matrices. A diagonal entry of 
Ky could be conductance; in that case the respective entry 
of Kj is unity. Alternatively, a diagonal entry of Kj could 
be a branch resistance; in that case the respective entry of 
Ky is unity. The components of the fault vector f G describe 
the faults in the branches. 

We consider three types of faults. First, a fault of switching 
element k, which is thought open but might be actually 
closed, is expressed in the form (6) as 

- jk = fc,k , (7) 

where the fault parameter f G ,k is a current through the 
switching element. The current should be zero if the switch- 
ing element is open. 

The fault of a switching element k, which is thought closed 
can be expressed in the form (6) as 

Vk = fc,k , (8) 

where the fault parameter f G ,k is a voltage drop at switching 
element. The voltage drop should be zero if the switching 
element is closed. 

Finally, consider the fault of the DC load. We assume that 
the load resistance is known. The fault makes load current 
in branch k to differ from what is given by BCE model (6) 

gkVk - jk = fo,k, (9) 

where g\ ; is the branch conductance and fc.k is a load current 
deviation from the model, which is indicative of the fault. 
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B. Observation equations 

The STA equations need to be complemented by obser- 
vation equations. We model sensor faults through additive 
offsets caused by these faults. 

We assume that the observations include the currents 
jmeas £ ft Ml measured by current sensors and voltages 
tt-meas £ ft Mv measured by voltage sensors. The respective 
observation equations have the form 

6meas — SyC T fv , (10) 

jmeas = S ij + fj , (H) 


(11), (12), and (6) and adding noises £ and rj we obtain the 
linear model of the form (1), (2). In this model x, y, A, and 
C are given by (14), (16); the faults are defined by (15) and 
the fault-related matrices B, D in (1), (2) have the form (17). 

For zero noises £ = 0, rj = 0 the STA equations (1), (2) 
make a system of total M y + 2 N B + Ni n equations in the 
2Ng + Ny unknown components of vector x. In a special 
case of M y = Njy — Ni n (where Njy — Ni n is the number of 
the edge nodes), the system is square and a unique solution 
can be found. This paper considers the redundant observation 
case of M y > Ny — Ni n . 


where fy and fj are fault offsets, Sy and Sj are selection 
matrices (only some nodes and branches are instrumented). 

The known voltage sources e srce £ :R Ms and grounded 
nodes e grn d £ 3? Mg (we assume that e gr d = 0) are counted 
among the observations. Though no on-line voltage sensing 
might be available for the source voltages, there might be 
an off-line knowledge of these. The respective observation 
equations have the form. 

tt-srce — Sge T fS: (12) 

tt-grnd. — Sg&i (13) 

where fg is the fault offset, Sg and Sc: are selection 
matrices. We do not consider ground voltage faults though 
it would be easy to add these to the model. 

The voltage offset fy, the current offset //, and the 
source (battery voltage) offset fg are unknown vectors of 
appropriate dimensions. 


C. STA model with faults 

To integrate the STA equations, introduce a state vector 

X £ Sft 2,V n+V,v alK j a ma trix A £ Sft 2N B+Ni„,2NB+NN as 
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Let / be a vector of all fault parameters. 

/= [fa fv fl (15) 


where Np = Np + My - 1- Mi + Mg. Introduce the 
observation vector y £ i R My , where M y = Mi + My + 
Mg + Mq, and the observation matrix C £ sft A L> 2 Wi+.'v,v ) 
The combined observations can be expressed as 
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Consider matrix B £ ^' 2N ri+N,, n ,.N F d escr ibing the fault 
impact on the states x and matrix D £ 9t N v ,N - F describing 
the impact of the sensor faults on the observation vector y 
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The faults enter equations linearly. After pulling together 
state equations (4), (5), (6) and observation equations (10), 


IV. Verification 

This section presents an application to ADAPT. We de- 
scribe a model of the testbed, then discuss simulation results. 

A. Modeling 

The ADAPT circuit is shown in Figure 1, see [23], The 
battery voltages V\ and If are assumed to be known. There 
are six branch resistances (loads) that are assumed to be 
known as well: these include two internal resistances of the 
batteries RBI, RB2; two DC load resistances RDC1, RDC2; 
and two equivalent resistances RAC1, RAC2, of the AC load 
branches (including the AC/DC inverters). The resistances 
of the circuit breakers and relays are assumed to be zero in 
the closed state and infinity in the open state. The circuit 
parameters were identified from experimental data. 



Fig. 1 . Data tags for the ADAPT cicuit 


The six current sensors in the circuit yield components of 
the current measurement vector j meas (11). There ten voltage 
sensors yield components of the voltage measurement vector 

&meas in (10). 

The branch constitutive equations (6) are defined by the 
states of the switching elements (six relays and six circuit 
breakers). These open/closed states are collected from the 
embedded switching element sensors in ADAPT. 

B. Simulation 

The incidence matrix of the circuit was defined in accor- 
dance with Figure 1. The source voltages and resistances 
used in the simulation are summarized in Table I. 
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Vi 

v 2 

RBI 

RB2 

RAC1 

RDC 1 

RAC2 

RDC2 

25.84 

24.83 

0.1 

0.1 

6 

10 

4 

20 


1SBEE I 

ADAPT SIMULATION PARAMETERS 


ESH141A, ESH144A 

ESH160A, ESH241A, ESH244A. ESH260A 

open 

closed 


TABLE 11 

Basic relay configuration in the simulation 


The simulations used the STA model (1), (2), (14), (16), 
(15), (17) with zero noises £ and //. In simulation, current and 
voltage measurements were ignored by assuming Mj = 0, 
My = 0, and eliminating the vectors j m eas > e meas , from 
(11), (10). The resulting system of equations has 2 Nb + 
Nn = 54 states x and 2 Nb + Ni n + Ms + Mq = 54 
equations, where Nb = 18, Nn = 18, Ni n = 12 , Ms = 2 , 
and Mq = 4. The Ms source voltages and Mq ground 
voltages (6 at all) define N^ — Ni n = 6 boundary conditions. 
For given switching element states, the resulting square 
system was solved to determine the state vector x. 

The switching element states define the BCE matrices 
K [ and Ky in ( 6 ), (14). The basic switching element 
configuration in the simulation is shown in Table II. All 
breakers are assumed closed. The relays are configured to 
connect Battery 2 to all DC and AC loads; Battery 1 is 
disconnected from the loads. 

The simulation input is the fault state /, which we 
generated as a sparse random vector. The switching element 
faults were applied by inverting the respective open/closed 
states before calculating the BCE matrices AT/, Ky , and the 
STA matrix A (14) in accordance with (7), ( 8 ). 

The load faults were applied by modifying the respective 
load resistances in BCE ( 6 ). The resistances shown in Table I 
were modified by a given percentage in the range from - 
50% (-100% corresponds to a short circuit) to +50%. The 
source faults were applied by modifying e src in (16). The 
source voltages shown in Table I were modified by a given 
percentage in the range from -50% (-100% corresponds to 
zero battery voltage) to zero (no change in the voltage). 

The current and voltage sensor measurements at the simu- 
lation output were modified by adding offsets proportional to 
the respective fault magnitudes. The fault offsets for current 
sensors were distributed in the range from -1A to 1A. The 
fault offsets for the voltage sensors were distributed in the 
range from -12V to 12V. 

In the simulations, all parametric faults were constrained 
to be at least 20 % of the respective maximum magnitude. 
This avoids small faults that are below the noise level. 

C. Diagnostics of simulated data 

The diagnostics algorithm uses model (1), (2). The main 
differences with the simulation model are as follows. 

• The diagnostics model includes faults. The fault model 
is described by (15), (17). 


• In the diagnostics model, the switching elements are 
assumed to be in the nominal state. If a fault is present, 
the actual state differs from the nominal. 

• Observations (2), (16) in the diagnostic model are the 
currents and the voltages distorted by noise 77 . 

• The observation and process noises £ and 77 in the 
diagnostic model ( 1 ), ( 2 ) are described through the 
inverse covariance matrices R and (). 

The dimension of fault vector f (15) is Np = Nb + 
My + Mi + Ms = 38. The fault estimates were computed 
by solving QP problem (3) and then thresholding the absolute 
values of the fault vector /. 

The problem (3) is transformed into a standard form of a 
QP problem with 130 decision variables. The solution using 
Mosek, see [20], takes about 20 msec on a 2 Ghz Wintel 
laptop computer. This is suitable for real-time implementa- 
tion and much faster than 0.5 sec sampling time for ADAPT 
data collection. 

In the simulation experiments, a series of fault patterns 
were seeded (input into the simulation) to generate the sensor 
data. The estimation algorithms are then applied and the 
diagnosis of the fault state is compared to the seeded faults. 
In great majority of the simulation runs, the estimation 
algorithms diagnose the seeded faults exactly. We will focus 
on the cases when the diagnosis was inexact. 

The first series of tests was to seed all faults, one fault 
at a time. There were 5 imperfect diagnoses encountered 
among the 38 total cases, one case for each fault component. 
While imperfect, these diagnoses were not, strictly speaking, 
incorrect. The imperfect diagnoses are discussed in Table III. 

In the second series of tests two faults were seeded at a 
time. About 200 two-fault cases were randomly generated out 
of the entire set of the 38-37/2 = 703 double fault cases. In 
a majority of cases, a perfect diagnosis was achieved. Each 
time one of the five faults from Table III was encountered, 
the diagnosis was imperfect as described in the table. In 
addition to that, combinations of two faults introduce new 
ways for the diagnostic algorithm to be in error that are 
listed in Table IV. While imperfect, most of these diagnoses 
are not incorrect. 

1) Circuit breaker of a disconnected battery fails open. The currents 
do not change, but disconnected voltage sensor now floats. This is 
misdiagnosed for voltage sensor fault. 

2) Load circuit breaker fails open. This is misdiagnosed for an offset of 
the load current sensor, which would yield the same measurements. 
The faulty part of the circuit is pointed out correctly. 

3) Same as 2, but for a different load branch. 

4) Fail closed relay is diagnosed correctly. In addition to that, two more 
spurious faults are reported. The observed change in the current 
balance could be attributed to each of the three reported faults. 

5) Battery output voltage is offset. The fault is detected correctly, but 
there is a false positive. The battery circuit breaker is misdiagnosed 
failing open. No current flows through the CB and it failing open 
cannot be excluded based on the data observed. See # 1 for more 
explanation. 

TABLE III 

Single fault diagnostics errors 
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1) A failed circuit breaker or relay are diagnosed correctly. A second 
fault in the circuit part unpowered by the first fault remains unde- 
tected. 

2) A fail closed relay creates a parallel connection of the two batteries. 
This fault is diagnosed correctly. A second relay failure is misdi- 
agnosed because the parallel connection of the batteries makes the 
currents through them very sensitive to the (close) battery voltages 
and (small) internal resistances. 

3) Resistances of the two DC loads are offset, but only one of the 
two faults is recongnized. The second fault is below the threshold. 
Decreasing the l\ regularization penalty would allow detecting the 
second fault but lead to false positives in other cases. 

TABLE IV 

Two-fault diagnostics errors 


In summary, the proposed approach works quite well in 
simulation. In majority of the cases the diagnosis is perfect. 
In great majority of all cases, the problematic circuit branch 
is diagnosed correctly. The remaining imperfect diagnoses 
make for less than 5% of all cases and are caused by high 
or low problem sensitivities to some of the fault parameters. 

The ADAPT circuit includes DC/AC inverters and AC 
loads. The inverter together with the AC load circuit can 
be modeled through an equivalent DC resistance. Such 
linearized model of inverter can be obtained by fitting experi- 
mental data on inverter draw current and supply voltage. The 
suitability of such a simplified inverter model is confirmed 
by the experimental results. 

V. Validation 

The ADAPT testbed is an instrumented EPS with two 
sets of computer controls. The data from the first set are 
called “Observer” data and are used for testing diagnostic 
algorithms. These data were used in this work. The second, 
independent set of controls, known as “Antagonist” is em- 
ployed in creating fault conditions for the Observer. These 
controlled fault conditions are known and used as a reference 
when evaluating diagnostic method performance. More detail 
on ADAPT can be found in [23]. 

Once verified in the simulation, the algorithms were inte- 
grated with ADAPT and validated in experiments. 

One data set used for the algorithm testing is displayed 
in Figure 2. The data was collected at 2 Hz rate (0.5 sec 
interval). In the experiment, all loads are initially unpowered. 
Then relays EY244A and EY260A close to connect Battery 
2 and to power AC Load 2. After a transient, the DC/AC 
Inverter 2 powers the AC load nominally. After several sec- 
onds, circuit breaker ISH262 fails (opens). This disconnects 
the load. The faulty circuit breaker sensor continues to show 
that it is closed. The described scenario is complemented 
by a voltage sensor fault (offset ramp). Figure 2 shows 9 
plots of selected observed signals, those with a transient and 
those later found faulty. The remaining 27 channels out of 
the 36 monitored show constant values, e.g., battery voltages 
or circuit breaker states (all closed) and are not plotted. 

Though the data contains transients, the diagnostic al- 
gorithm assume a steady state model at each step. The 
steady-state model used for diagnosis neglects the transients. 
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Fig. 2. Observed experimental data 
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Fig. 3. Estimated fault states for the experimental data 


Despite that, the algorithms are able to detect and diagnose 
the seeded faults. The results are illustrated in Figure 3, 
which shows the fault estimates corresponding to the signals 
in Figure 2. The diagnosis is perfect with two exceptions. 
First, the inverter warm-up transient is diagnosed as a fault. 
Second, it misses the sensor fault in the circuit branch that 
becomes unpowered. 

The inverter model assumes a steady-state condition. How- 
ever, at the start of the inverter transient at 27.5 seconds, the 
inverter takes approximately 6 seconds before power appears 
at the output. Then, it takes an additional 14 seconds before 
all the loads are turned on and the equivalent resistance 
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Time: Event Description 

t=24: Control action Relays EY244A and EY260A are closed to 
power the load 

t=27.5: Start of INV2 INV2 Input CD ISH262 is suspected open; 
transient INV2 Load is suspected to be off-nominal 

t=34: Sensor fault oc- E242 voltage offset starts ramping down from 

curs zero 

t=46.5: Sensor fault de- E242 voltage offset is detected during the 
tected inverter transient 

t=47.5: Nominal regime Inverter reaches a steady state 

t=59.5: Circuit breaker Circuit breaker ISH262 opens and cuts off 

fault occurs the branch with the DC/ AC inverter load and 

E242 voltage sensor 

t=60: Circuit breaker Circuit breaker ISH262 failure is detected; 

fault detected, sensor E242 voltage offset fault in the cut-off branch 

offset detected stops being detected 

t=65.5: Relay fault de- The large offset of voltage measurement E242 

tected, sensor fault un- is now interpreted as a voltage across relay 

detected, false detection EY244A in the disabled branch indicating 
of a second relay fault that the relay failed open. 

TABLE V 

Timeline of events in experimental data and fault estimation 


reaches the full RAC2 value that is used in the model. 

The lower right plot shows an estimate of the fault 
current (deviation from the nominal current) for the DC/AC 
Inverter 2 with the AC load. Switching of the relays and 
inverter excites very large current transients. These causes 
false positives for current sensor faults. The middle plots in 
Figure 3 show three false positive occurrences for one sensor 
and two for another, one time-sample duration each. The 
timeline of the events and the diagnostic estimation results 
with a commentary are detailed in Table V. 

VI. Conclusions 

The proposed diagnostics approach works very well for 
the DC power system applications exemplified by ADAPT. 
The diagnosis is accurate when this can be expected. The 
algorithm has a few well-defined parameters that need to 
be set up and takes milliseconds to compute results. The 
approach has been demonstrated to work well with exper- 
imental data despite ignoring the transients and modeling 
a complex circuit of DC/AC inverter with a load by a 
single DC resistor. The approach is clearly suitable for on- 
line implementation in practical monitoring applications for 
electric power systems. 
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